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ABSTRACT 


This thesis supplies a statistical and economic tool for 
analysis of the failure characteristics of one typical piece 
of equipment under evaluation: a cam-driven reciprocating pump 
used in the submarine's distillation system. Comprehensive 
statistical techniques and parametric modeling are employed to 
identify and quantify pump failure characteristics. Specific 
areas of attention include: the derivation of an optimal 
maximum replacement interval based on costs, an evaluation of 
the mission reliability for the pump as a function of pump 
age, and a calculation of the expected times between failures. 
The purpose of this analysis is to evaluate current 
maintenance practices of time-based replacement and examine 
the consequences of different replacement intervals in terms 
of costs and mission reliability. Tradeoffs exist between 


cost savings and system reliability that must be fully 
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understood prior to making any policy decisions. 
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EXECUTIVE SUMMARY 


Today's Submarine Force is faced with many new 
challenges. These include an ever-expanding role in Joint and 
Maritime operations, despite a diminishing operating budget 
and a shrinking force level. Submarine maintenance is one 
key area that must be carefully evaluated to ensure maximum 
operational readiness while working within the confines of 
budgetary constraints. A Maintenance Effectiveness Review, 
MER, is currently in progress to evaluate the applicability 
and effectiveness of all aspects of submarine maintenance. 
Any comprehensive evaluation of maintenance practices, 
however, requires a detailed analysis of equipment performance 
in terms of failure characteristics and life expectancy. 

This thesis supplies a statistical and economic tool for 
analysis of the failure characteristics of one typical piece 
of equipment under evaluation: a cam-driven reciprocating pump 
used in the submarine's distillation system. The purpose of 
this analysis is to evaluate current maintenance practices of 
time-based replacement and examine the consequences of 
different replacement intervals in terms of costs and mission 
reliability. Tradeoffs exist between cost savings and system 
reliability that must be fully understood prior to making any 
policy decisions. 

The analysis examined the failure characteristics of a 
sample set of 61 pumps installed between May 1987 and December 
1993. The data consisted of 140 failures, of which 14 


required pump replacements, and 2020 total months of pump 
operation. The primary source of the data is the Navy's 3-M 
system. Pump failures are classified as repairable or non- 
repairable and thus require replacement. The probability that 
a given pump failure is repairable was evaluated to be 0.9. 
Additionally, a new pump will experience on average 10 
failures cycles prior to being replaced. 


Trending analysis indicated an increasing failure rate 


with pump age; a sign of system wearout. The increasing 
failure rate, likewise, indicated that a time-based 
replacement maintenance policy may be warranted. A 


Nonhomogeneous Poisson Process (NHPP) was chosen to model the 
pump's failure characteristics. Maximum likelihood was then 
used to estimate the NHPP's parameters. 

The use of a stochastic modelling allowed a comprehensive 
evaluation of the current maintenance policy of time-based 
replacement of the pump at a periodicity of 36 months for 
Trident submarines and 60 months for non-Trident submarines. 
The pump is the same for both platforms. The optimal 
replacement interval based entirely on minimizing the long-run 
average system costs was determined to be 111 months for an 
average lifecycle cost of $2226 per month. This periodicity, 
however, resulted in undesirable mission reliability. The 
resulting probability of the pump completing a mission of 
three months just prior to replacement was only 0.23 with 1.47 
expected pump failures during the mission. A replacement 
periodicity of 36 months resulted in a mission reliability of 
0.76 with 0.27 expected failures at an average lifecycle cost 


of $3674 per month. The 60 month periodicity hada mission 


reliability of 0.63 with 0.46 expected failures at an average 
cost of $2611 per month. 

The results of this analysis indicate that a 60 month 
replacement interval may be acceptable based upon reliability 
reliability requirements. The periodicity of 36 months, 
however, may be excessive and warrant extension to 60 months. 
This thesis does not attempt to quantify any minimum 
reliability requirements. It is evident, though, that a 
replacement schedule based strictly on economic considerations 
is unsatisfactory from a reliability standpoint. This 
illustrates on of the most powerful uses of stochastic 
modelling: the ability to predict and evaluate the 
consequences of maintenance policy decisions on system 
performance. 

The scope of this thesis has implications beyond the 
operation of this one pump. This case study of the cam-driven 
reciprocating pump illustrates the type of analytical 
techniques necessary to perform a comprehensive evaluation of 
a shipboard system's performance. The goal is to provide 
decision-makers with the in-depth statistical foundation to 
make sound decisions on maintenance policy; decisions that 
directly affect the readiness and ability of the U.S. 


Submarine Force to assume its expanding mission. 


INTRODUCTION 


A. BACKGROUND 

The end of the Cold War and growing concerns over the 
national debt have forced the Navy, and specifically the 
Submarine Force to reevaluate its current operational policy. 
One of the areas being closely evaluated is that of Submarine 
Maintenance. RADM (sel) R.E. Frick, NAVSEA Deputy Commander 
(Submarines Directorate), recently directed an evaluation of 
all submarine maintenance performed at Depot (shipyard), 
Intermediate Maintenance Activity (IMA), and Ship's Force 
levels. This Maintenance Effectiveness Review (MER) is a 
continuous process to evaluate the maintenance needs of the 
Submarine Force while working within the fiscal constraints of 
a shrinking budget. NAVSEA PMS390, Submarine Monitoring, 
Maintenance, and Support Program Office (SMMSO) is the lead 
activity in charge of coordinating this review. 

The review process requires an in-depth review of the 
applicability and effectiveness of all maintenance actions 
performed by the submarine force. The goal is to eliminate 


redundant and ineffective maintenance while maintaining 


optimal materiel readiness, within a limited budget. The 
concepts of Reliability-Centered Maintenance are being used to 
evaluate the time-based and "fix when fail" policies 
historically used by the Navy.[Ref.1] The key to a successful 
shift in strategies is the ability to monitor equipment 
conditions and make reliable maintenance decisions based on 
real-time data and predictive analysis. 

The goal of this thesis is to address the analytical needs 
of the Submarine Force necessary for a thorough assessment of 


current maintenance policies. 


B. CURRENT ANALYSIS 

The SMMSO organization consists of a staff of approximately 
one hundred personnel based in Washington, D.C., as well as 
monitoring teams stationed at each submarine base. The focus 
of SMMSO's work is submarine equipment performance monitoring 
and maintenance recommendations. On a periodic basis, SMMSO 
site teams perform various performance tests on a submarine's 
systems. These tests involve everything from the vibration 
monitoring of a pump to oil analysis of shaft lubricating oil. 
The results of these tests, as well as any significant 
information on past problems, are analyzed by a specific 


engineer who is SMMSO's system expert on that equipment. This 


information is used to predict future equipment performance, 
make recommendations on corrective actions, and make 
maintenance deferment and equipment replacement decisions. 
In the past, the scope of the statistical analysis was 
focused on short-term predictions, i.e., refit to refit, which 
was at the time, adequate to formulate maintenance policy 
given the large Defense Budget and high priority on submarine 
readiness. These policies have been very conservative in 
nature and possibly more restrictive than necessary. In view 
of the shrinking budget, these policies should be reviewed to 


ensure the optimal use of availible resources. 


C. PROBLEM STATEMENT 

The specific area of analysis to be addressed in this 
thesis is the performance, from the standpoint of reliability, 
availability, maintenance, and replacement needs of a 
particular cam-driven reciprocating pump. The Class 
Maintenance Plan (CMP) of this pump is currently being 
evaluated for effectiveness and applicability. Mr Richard 
Youngk is the SMMSO systems engineer conducting a performance 
analysis of the cam-driven reciprocating pump. This thesis is 
in conjunction with his efforts and builds upon his research 


and analysis. [Ref .2] 


The selection of this particular pump for the study is in 
itself not significant, but the pump is typical of the systems 
being analyzed. The analysis presented in this thesis focuses 
ina broader context on the specific questions being asked 
more broadly. This thesis proposes a process or method to 
systematically analyze the pump data presently available. It 
is believed that this procedure can then be adapted for use to 
analyze other systems. 

Specific questions to be addressed in the analysis are: 

1. What is the expected number of pump repairs and 
replacements required over a certain time interval? 


What is the associated expected repair cost and 
downtime? 


2. Does the failure rate change over the life of the pump? 
3. What is the optimal equipment replacement interval? 


4. Can a useful stochastic model of pump performance be 
constructed? 


D. FAILURE 
Before proceeding further with the analysis, it is 


necessary to define what is meant by the term: failure. 


A failure is classified as an event or inoperable state, 
in which any item or part of an item does not, or would not, 
perform as previously specified. [Ref.3] This paper further 
classifies failures as being either repairable or non- 
repairable. A repairable failure implies that the system can 
be returned to operating specification with the replacement of 
only a small fraction of the system's parts and in a short 
period of time. Non-repairable implies the pump must be 
completely replaced, or else an extensive overhaul is 
necessary to restore pump operation within allowable limits. 
In the context of this paper, pump replacements are limited to 
replacements resulting from a non-repairable failure and not 


replacements based upon any other criteria. 


E. FAILURE MODELS 

Two basic models often used in evaluating repairable 
systems are the Homogeneous Poisson Process (HPP), and the 
Nonhomogeneous Poisson Process (NHPP). Both the HPP and NHPP 
are counting processes used to model the number of component 
or equipment failures occurring over the life of the 
equipment. Figure 1.1 shows the life-cycle of a piece of 
equipment, with t, representing the time to the it? failure, 


and x; representing the time interval between the (i-1)* and 


it’ failures, i.e. xX, = ty-ty.,. To be more explicit, > a 
represents the time from the (i-1)* to i‘* failure for item 
(here pump) k. Replacement of pump k means that the 
repairable-failure generating Poisson process starts again 


from scratch. 


Figure 1.1: Model of Pump Life 


Now let N(t) equal the number of failures occurring prior 
to pump age t. Then N(t) is considered a counting process, 
{N(t), t20}, if the following conditions hold: 


N(E) 20. 

N(t) is an integer. 

If s < t, then N(s) < N(t). 

For s < t, N(t) - N(s) represents the number of events 


occurring in the interval (s,t). 


BWh PH 


1. Homogenous Poisson Process 
The Homogenous Poisson Process is a special counting 
process which is commonly used to describe systems with a 
constant failure rate A. More precisely, the counting 
process, {N(t), t20}, is said to be a Homogenous Poisson 
Process with rate A, A > 0 if: 

1. N(t)=0. 

2. The process has independent increments, i.e., the 
number of events occurring in disjoint time intervals 
are independent. 

3. The number of events in any interval of length t=tT;-T, 
is Poisson distributed with mean At 

i.e. for n = 0,1,2,..., 


(A(T,-T,)) * 
n! (1.1) 


Pr[N(t,)-N(T,)- n] - eh ltyrty) 


Ti<T yr n-O,1,... 


Furthermore, the conditional probability that the system will 
survive to time 1, given that it is operating at time 1, 


denoted by R(t;,1T;), 1s 


-A(t,-t,) 
R(t ,t Jee ae 


(1.2) 


TST, 


2. Nonhomogeneous Poisson Process 
The Nonhomogeneous Poisson Process, also called a 
nonstationary Poisson process, is a counting process similar 
to the HPP except that the rate function is not constant, but 
a function of system age, t, denoted by A(t). The x,'s, i.e. 
the times between equipment failure, are not necessarily 
identically distributed. This mathematical structure can 
represent an increasing failure rate over time resulting from 
equipment wear. Pump data indicate such behavior. 
Specifically, the counting process, {N(t), t20}, is said 
to be a Nonhomogeneous Poisson Process with age-dependent 
repairable failure rate A(t) if: 
1. N(t) = 0. 
2. N(t) has independent increments. 
3. The number of events, i.e. pump failures, in the 


interval (1T;,1T;) is Poisson distributed with 
mean m(1%;,T3), 


y 
m(% pt )= [r@)as. (1.3) 


YW 


The probability of n failures in the time interval 


(Ti, 7,5) is 


: . - ~on(4,,5)) (m(«,,t,))” 
PrIN(s,)-N(x)- n] = — (1.4) 


T<Tp n-0,],... 


The reliability of the system at any time, t, depends on 
the age at which the most recent failure occurred. Suppose 
that a failure actually occurs at time t = 1;, where 1; is the 
age at which the jt failure of the incumbent system (pump) 
occurs. Then the probability that the system will not fail 
for at least age T, + x units of time is 


Pr[t zx|t pt] - e ntte) (1.5) 


pr'y 


II. MODEL DEVELOPMENT 


A. DATA 

The data used in this study consist of a sample set of 
sixty-one pumps and 140 failures over the observation period 
beginning in May 1987 and ending in December 1993. The final 
observation was either a failure and replacement prior to 
December 1993, or the pump was last observed still in 
operation. This sample set does not include all pumps in 
operation, but a subset for which there exists reliable dates 
for pump installation. Additionally, the pumps were carefully 
screened to ensure that all pump replacements resulted from 
non-repairable failures. The installation and failure times 
are rounded to the nearest month and henceforward all 
references to time and age will be in months. The time in 
service is adjusted to remove any inactivation period of two 
or more months in length from the pump's operating age. The 
pump itself has no runtime meter. Engineering logs monitor a 
pump's operation and would provide a more accurate measure of 
the actual operating time for each pump. However, it is not 


feasible to accumulate and analyze this information ina 
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reasonable time frame. Failure and replacement data originate 
from the Maintenance and Material Management (3-M) system 
which reports routine maintenance via OPNAV Form 4790/K, 
Casualty Reports (CASREP), the Submarine Maintenance, 
Engineering Planning and Procurement (SUBMEPP) system, and 
input directly from the Fleet. 

Mr. Richard Youngk conducted a Failure Modes and Effects 
Analysis, FMEA, ‘to identify the critical failure modes 
characteristic to the pump and characterize the required 
repair actions. Information on the types of failures, their 
frequency, the repairs requirements, and the trends for the 
different failure modes is essential elements in any analysis 
of equipment performance. This thesis focuses specifically on 
the occurance of failures and the required repair action. 
Failures are not distinguished, but treated as one entity. 
The lifecycle cost in terms of material and labor is also 


examined. Appendix A is a copy of the actual data analyzed. 


B. ASSUMPTIONS 
Several assumptions were made regarding the pump failure 
data. These assumptions are: 


1. The submarines in the sample set have relatively 
similar operating cycles, i.e., pumps on different 
submarines will face the same operating conditions 
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and undergo roughly the same number of hours of 
operation for the same time period. 


2. Every component failure causes equipment failure. 


3. Failures are immediately evident. 
4. Pumps are only repaired or replaced at failure and not 
in anticipation of failure. In actuality this is not 


always the case, but the data set was cleansed to ensure 
only replacements related to failures are counted. 


5. Consecutive failures on an individual pump are 
independent. 


6. A repair returns the pump to full operation, but does not 
necessarily restore it to a "good as new" condition. 


7. Equipment repairs consume no appreciable time. 
8. A replacement constitutes the installation of a new 
pump or a complete overhaul of the current pump. 
The accuracy of these assumption and the accuracy of the data 


will be discussed further in Chapter IV. 


C. ANALYZING THE DATA FOR TIME DEPENDENCY 

One of the initial steps in the analysis is to determine 
if the data are compatible with an increasing failure rate 
over time, i.e., with aging or wearout. Evaluating the time 
intervals between failures provides some insight into this 
matter. For a constant failure rate over time, the intervals 


between successive failures should be very similar; with all 
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interval times following an exponential distribution with the 
same mean. Care must be taken, however, to ensure all data 
has been taken into account. Invalid conclusions can be 
drawn, for example, by comparing the mean time to first 
failure with the mean interval time between the first and 
second failure if all pumps have not failed at least twice. 
Here the existing life of the pumps with only one failure is 
not accounted for in the mean for the second interval. 

A better approach is to compare the mean interval time 
between successive failures for pumps with like failure 
numbers. Table 2.1 displays the mean and standard deviation 
(STD) for the interval times between failures for pumps with 
at least two failures, with at least three failures, and so 
on. In all cases the overall trend is a decreasing interval 
between failures over time, indicating an increasing failure 
rate. The correlation coefficient between successive failure 
times also provides evidence of an increasing failure rate. 
Table 2.1 also shows the correlation between the interarrival 
times of successive failures. The correlation between X, and 
X,, the interval times for failures one and two, and between 
X, and X, are not significant. The negative correlation 
between X, and X, and between X, and X,, however, indicates a 


decreasing interval time. 


TABLE 2.1: COMPARISON OF INTERVAL TIMES BETWEEN FAILURES 


A more rigorous statistical test for distinguishing between 


a constant failure rate, i.e., a HPP, anda monotonic trend is 
the Laplace Test. The Laplace Test is based upon the HPP 
property that given n arrivals in time (0,F), the unordered 
times of arrival, denoted by T,,T,,...,T,, are the ordered 


Statistics from an independent uniform random variable on 


(0,F). 
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Therefore, the test statistic for the Laplace Test, 


has approximately a standard normal distribution. 
(1955) 


significance level for n 2 4. [Ref.4] 


(2.1) 


Bates 


showed this approximation was adequate at the 5% 


A slight modification is required if the system is observed 


until a specific number of failures occur. In this reference 


case F = T, and the modification is 


(2.2) 
1 
12(n-1) 
Cox and Lewis (1966) showed that the Laplace test is 
optimal for the NHPP with rate 
A(t=e «PF 
(2.3) 


-0o <a, BP <~,t2 0 


in testing the hypothesis of a constant rate, 1,@.; Hot B=0; 
against the alternate hypothesis of a trend, H,:B#0. 

The test can be further modified to examine a series of k 
independent systems with the same value of B. The new test 


statistic is, [Ref.5] 


gk On k 
Z 1 


DBD DE Fad nf, 
y . Let 27 (2.4) 


1 gE : 5 
[a5 vf 


This procedure was applied to 13 pumps meeting the criteria 
of four or more failures. Table 2.2 gives the sample set used 
for the test. The evaluation of two pumps, number 57 and 60, 
required the modification of Equation 2.2 since the 
observation period ended with a failure for both 
pumps. 

The resulting value of the test statistic, U, was 
U =2.016 which resulted in a p-value of .04 for the two-tailed 
hypothesis test. Thus, the null hypothesis can be rejected in 
favor of a changing failure rate at the conventional 5% 
significance level. Further, the positive value of the test 
Statistic indicates an increasing failure rate over time. 
This test included only 13 pumps of the original data set of 


61 pumps. It is reasonable to assume, however, that this 
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TABLE 2.2: LAPLACE TEST FOR DATA TREND 


LAPLACE TEST FOR A TREND IN THE FAILURE RATE 
H,: B = 0, Hy: B # 0 


NUMBER 
OF 
FAILURES 


SUM OF 
FAILURE 


FINAL 
OBSERVATIO 
N 


STATUS 
C=RUNNING 
R=REPLACED 


TEST STATISTIC VALUE 
P- VALUE 


NOTES: * THE TOTAL NUMBER OF FAILURES IS MODIFIED FROM n; 
TO n;, FOR PUMPS 57 AND 60 SINCE THE FINAL 
OBSERVATION WAS A REPLACEMENT. 
** PUMP 50 IS NOT USED SINCE THE ADJUSTED NUMBER OF 


FAILURES IS 3. 


result is representative of the entire sample set. The 
Laplace Test was in fact performed on all pumps with at least 
one failure. This resulted in a value U = 1.912 and P = .06. 
Although possibly not as accurate as for the case with ne: 4, 
it does seem to indicate the overall trend is an increasing 


failure rate. 


D. CHOOSING A MODEL 
Given that the rate of failure occurance seems to increase 
with pump age, the next step involves determining whether a 
mathematical model can be constructed to accurately simulate 
the occurrences of failures. The nonstationary trend of the 
data indicated that the times between successive repairable, 
i.e., non-fatal, failures were not identically distributed. 
A particular Non-homogeneous Poisson Process (NHPP) has been 
selected to model the occurrences of failure over the life of 
the pump. This approach was chosen based on the fact that 
the NHPP assumption, besides being mathematically 
tractable, is a good representation of many systems in the 
real world because it is a consequence of the "minimal 
repair" assumption. "Minimal repair" implies that the 


repair involves the replacement of only a small fraction 
of a system's constituent parts. [Ref.6] 


18 


The following parametric failure rate, A(t), was used to 


model the increasing trend in failure over time: 


A(t-e © 
(2.5) 


-0 <a, BP <~,t2 0 


This form of failure rate was introduced by Cox and Lewis 
for the analysis of failure data for aircraft air-conditioning 
equipment in 1966.[Ref.7] Ascher and Feingold (1966) used 
this model to analyze the performance of submarine main 
propulsion diesel engines. This model has the advantage that 
since equation 2.5 is positive for all value of o and 6, no 
nonlinear restrictions are necessary for the estimators of 
these parameters. [Ref .8] 

Unlike the Cox and Lewis situation, however, each failure 
may result in either a repair or a replacement. Therefore 
each failure is subject to a probability of being repairable, 
denoted by p(t,,), and a corresponding probability of being 
fatal and requiring pump replacement, denoted by q(t,;,;) with 
q(tz,5) = 1 - p(ti,;). In fact, some repairs may be relatively 
easy while the platform is on a mission; others have greater 
operational impact in that they may require mission 


termination for repair to be made. 
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E. MAXIMUM LIKELIHOOD ESTIMATION 

The parameters a and B have been estimated using the method 
of maximum likelihood. The maximum likelihood estimate, 
abbreviated as MLE, for 8 is found first, then a is found. 
A numerical solution is necessary; there is no simple closed 
form solution for 8. 


Before deriving the MLE, it is necessary to define the 


following: 
k = index for the number of pumps in the sample set 
j = index for the individual pumps, (= dwar 
i = index for the failure number for pump j, i = 1,2,..n; 
n; = the number of failures for pump j 
f. = time of the last observation for pump j 


*Note: time is equivalent to pump age in this analysis* 
t,, = time of failure i for pump j 
p(t,;) = probability that the failure occurring at time 
t;,; is repairable. 
probability that the failure occurring at time 
t;,; is not repairable and therefore requires 
pump replacement. q(ti,;)= 1=plty,5).. 
F,; = indicator variable indicating the status of the last 

observation 


q(ti,5) 


F; = 1 for non-repairable failure 


F; = 0 for pump last observed operating 
A(t,;;) = rate of occurrence of pump failures for 
pump age tj,,;- 
The failure may or may not be repairable. 
A(t,,;) = the integrated value of the failure rate from 
pump installation to age t;,;; the expected number 


of failures(repairable) up to age ¢t;,; 


8 = vector of parameters, 8 = (a, 8B) 
L(@) = Combined Likelihood Function for pump failures 
L.(@) = Conditional Likelihood Function for pump failures. 


Cc 


NOTE: If the last observation is a failure, then Bee agg 
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Appendix B contains the complete derivation of the maximum 
likelihood estimators for a and $ and the _ following 
explanation will only deal with the final results of that 
derivation. 

The combined likelihood function, L(8), for all pump 


failures reduces to: 


Say k k | gt.) 4; ko 
L@) e? (bef - HT«,). (2.6) 
jul $l j Fleet 


p(t, ,) 


Substituting the A(t) described by Equation 2.5 results in the 


following likelihood function: 


a 


a qt, )\7 en 
L@y A wie! 57 B | | LE 1-<,). (2.7) 


ja ny java 


The conditional distribution of the observations given n, 
events for pump j is formed by dividing the combined 
likelihood function by the marginal probability of pump j 
having n, total failures. Taking the logarithm of this 
conditional distribution produces the conditional log 


likelihood function. 


gH rk ia k 
log LOB Yt, 4 td, login,!) dog By n, -Y njog(e™4-1) (2.8) 
fl él 7) fl rl 
The maximizing value for B is then obtained by maximizing 
log [L,(6)]. This can be accomplished by setting the 


derivative of the Conditional Log Likelihood Function, with 


respect to B, equal to zero and solving for the root, oe 


k 
ae) (B)] ES, @. FLO 6. (2.9) 
jl tel fl B jel 1-e %y 


o~ 


Since no closed form solution for B— exists, a numerical 
method is necessary to solve for the root. First, d[log 
L.(6)]/dB is plotted over a range of B for a visual 
approximation of the equation's root, Figure 2.2. This "best 
guess" is then used with Newton's Method to find a more 
precise value for the root. The resulting maximizing value of 


Bis § = 0.02258. 
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MAXIMUM LIKILIHOOD ESTIMATION OF BETA 
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Figure 2.2: Visual Approximation of the MLE 
for 8B 


A similar process is used to find the value of a. The 
conditional likelihood function, however, cannot be used to 
solve for a as it was for B since the conditioning process 
removes @ aS a parameter. Therefore, L(6) must be used to 
find the maximizing value of a. Taking the logarithm of L(6), 
differentiating with respect to a, and setting the equation 


equal to zero results in the following expression: 


dflog LO) ys, 8s May. 
a ye ee (2.10) 


Solving the equation in terms of a results in the following 


expression: 
7 k 
By’ 2, 
ms = yl 
& - h}-— (2.11) 
¥ @ 1) 
rl 
Inserting the maximizing value of § = 0.02258 into 


Equation 2.11 produces the resulting maximizing value of a 
is @ = -3.189. Appendix C contains the MathCad 3.1 program 
of Newton's Method and the calculations used to determine the 
parameter values of § and é. 

The resulting maximum likelihood estimator for the 


failure rate of the pumps is: 


hina Oe tS 0, (2.12) 
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The confidence intervals for & and § are obtained by 
inverting the observed information matrix to form an estimate 
of the variance-covariance matrix for & and @. The 
unconditioned likelihood function, Equation 2.7, is used for 
calculating the confidence intervals since conditioning 
removes @ aS a parameter. The resulting 95% Confidence 
Intervals are given in Table 2.3. Note that zero is not 
contained in the interval for $8, giving further indication of 
an increasing failure rate. Appendix C contains the Mathcad 


3.1 program used for the confidence interval calculations. 


TABLE 2.3: CONFIDENCE INTERVAL FOR MLE OF & AND B 


PARAMETER 


F. EVALUATION OF MODEL FIT 

The use of any parametric model to quantify a system's 
behavior must also involve a goodness of fit analysis to 
determine the adequacy of the model. Two methods are used to 
assess the suitability of the afore mentioned NHPP as a model 


for the failure process. One method is graphical and the 
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other is a formal statistical procedure using Pearson's Chi- 
Squared Test. 

The first method compares the an empirical failure rate 
taken over equidistant intervals to the model's failure rate 
for the corresponding times. The interval length selected was 
eight months. During each interval the total pump exposure 
was calculated as the sum of the number of pumps operating 
during each month of the interval. Similarly, the total 
number of failures for each interval was determined. The 
interval failure rate is simply a point estimate formed by 
dividing total failures by total exposure. The interval 
midpoint was used to calculated the model's failure rate, 
A(t) <e"8t. Table 2.4 contains the calculated rates and Figure 
2.3 is the corresponding graph. 

The model closely approximates the interval failure rate 
up to the last two intervals, i.e., intervals 49-56 and 57-64 
months. This region indicates a decreasing failure rate. One 
should note, however, that the total exposure for these 
regions is relatively small compared to the other intervals. 
The small exposure, therefore, may not be representative of 
the overall trend. Further, one must recall the nature of the 
equipment being evaluated. The system is a cam-driven 


reciprocating pump; a mechanical system already noted to 
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TABLE 2.4: COMPARISON OF INTERVAL FAILURE AND THE 
MODEL FAILURE RATE 


MODEL 
FAILURE 


INTERVAL 
FAILURE 
TE 


OBSERVED 
FAILURES 


TOTAL 
EXPOSURE 
(MONTHS ) 


INTERVAL 
(MONTHS ) 


RA 
1 


posves [ze totam | 0.16098 | 


INSTANTANEOUS FAILURE RATE COMPARISON 


° 
an 
wo 


FAILURES PER MONTH 


OBSERVED --@-— NHpp MopEL 


Figure 2.3: Observed Interval Failure Rate vs. the 
Model Failure Rate, A(t)= e7 
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exhibit indications of wearout. Thus, the system is highly 
unlikely to experience any reliability improvement at this 
point in life. The system engineer at SMMSO independently 
reached a similar conclusion while doing a parallel study, 
i.e., that the observed decreasing failure rate late in life 
was a statistical fluxuation, resulting from sample size, and 
not truly typical of system performance. Of course it is 
possible that pumps on some vessels will actually improve with 
age, up to a point. This feature remains open for further 
investigation. 

The Pearson Chi-Squared Test is a formal statistical test 
for goodness of fit. Here the expected numbers of failures 
under the NHPP assumption are compared with the observed 
numbers of failures. The same eight month intervals are used 
with the exception of the last two which are combined to 
provide at least five failures per interval. Table 2.5 
contains the results. The chi-squared value is x? = 9.267 
with 7-3=4 degrees of freedom. This produces a p-value of 
0.05 <p < 0.1. Such a p-value is not ordinarily considered 
to indicate significant departure from the basic hypothesis, 
in this case the model. 

Note the large contribution of the last interval. The same 


arguments as above can be made for this deviation. Even so, 
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one cannot reject the possibility that the data follows the 
NHPP model at a 5% significance level. Therefore, the model's 
representation of the data, although not perfect, appears 
acceptable. Note further that acceptance of the model isa 


conservative step in predicting future failure characteristics 


for the pump. 


TABLE 2.5: CHI-SQUARED GOODNESS OF FIT TEST 


INTERVAL TOTAL OBSERVED EXPECTED CHI 
(MONTHS ) EXPOSURE FAILURES FAILURES SQUARED 
(MONTHS ) 


CHI SQUARED VALUE: x? = 9.267 with 7-31 
P-VALUE: .05 < p< .1 
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III. DEVELOPING A REPLACEMENT POLICY 


A. BENEFITS OF STOCHASTIC MODELING 

The benefit of developing a stochastic model of the cam- 
driven reciprocating pump's failure characteristics is that 
pump reliability is now described by a simple mathematical 
expression that describes likely future behavior. One might 
argue that empirical data would provide a_ better 
representation of actual pump performance. This is true for 
systems with an extensive data base on lifecycle performance. 
Such an extensive data base, however, does not exist for the 
cam-driven reciprocating pump. In fact, very little 
information exists for pumps over 60 months age. The lack of 
such data demonstrates the need for a probabilistic model to 
predict future behavior. The mathematical model allows easy 
calculation of the expected number of pump failures, and its 
availability for different pump ages if downtimes are also 
modeled. It also expedites determination of an optimal 


replacement interval based on costs. 
This chapter will discuss the merits of time-based 


replacement and evaluate the applicability and effectiveness 
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of a time-based replacement policy for the cam-driven 
reciprocating pump. The model and results developed in 


Chapter II provide the basis for the evaluation. 


B. CURRENT MAINTENANCE POLICY 

The cam-driven reciprocating pump is currently on a time- 
based replacement schedule. The schedule varies, however, 
between platform types, specifically between Trident and non- 
Trident submarines. The same pump is used in both platforms, 
although slight differences in the systems exist. Current 
replacement intervals for Trident and non-Trident assets are 
approximately 36 and 60 months, respectively. 

The time-based replacement policy involves the replacement 
of the cam-driven reciprocating pump at a predetermined age 
regardless of the current material condition of the pump. The 
value of this time-based replacement in maintenance planning 
is obvious. A planned maintenance evolution allows for the 
scheduling and prepositioning of parts, personnel, and support 
facilities to minimize system downtime and thus tends to 
minimize total system (submarine) nonavailability. It also 
precludes possible extensive and expensive failures that may 


occur at later ages. 
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This reasoning has led to the frequent and possibly 
excessive use of time-based maintenance by the Navy. 


The Navy uses time-based maintenance for components, 
equipment and systems ranging in complexity from oil 


filters to propulsion gas turbines. Most of the 
maintenance action in Class Maintenance Plans are based 
on engineering time based periodicities. RCM 


[Reliability Centered Maintenance] requires that these 
intervals be adjusted based on equipment performance and 


failure rate. Most time-based 
overhauls/refurbishments/replacements are also expensive. 
CBM [Condition Based Maintenance], applied where 


appropriate, will greatly reduce the number of time-base 
repairs and overhauls conducted. The key to successful 
implementation of CBM is application of the proper level 
of monitoring, evaluation and trending for each piece of 
equipment. [Ref .9] 


The challenge, therefore, is to evaluate the appropriateness 
and effectiveness of the current time-based maintenance 
practices and to optimize the interval used. The following 
issues must be addressed for the cam-driven reciprocating pump: 
_1. Is this pump an appropriate candidate for time- 

based replacement? 

2. If time-based replacement is warranted, what is 
the optimal replacement interval to minimize 
cost? 

C. CRITERIA FOR TIME-BASED REPLACEMENT 

The goal of time-based replacement is to improve both the 

current and long-run operating state of the system through 


preplanned maintenance actions. The meaning of "operating 
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state" is dependent upon the nature and mission of the system in 
question and upon the objectives of the policy makers. Thus, 
"operating state" refers to operating cost, system reliability, 
maintainability, and so on. 

MIL-STD-2173(AS) provides the following guidance on the 
applicability of time-based tasks, also referred to as hard time 
tasks. Here the time-based task is considered to be of two 
types, scheduled rework or scheduled discard. A reworking task 
is analyzed if reworking promises to restore the item to an 
acceptable level of failure resistance; otherwise, the discard 
task is analyzed. 

The applicability criteria for time-based (hard time) tasks 
are as follows: 


1. The item must be capable of having an acceptable 
level of failure resistance after being restored 
(for rework task). 


2. The item must exhibit wearout characteristics, which 
are identified by an increase in the 
conditional probability of failure with increasing 
usage(age). This property can lead to establishment 
of a wearout age (for rework tasks) or a life-limit 
(for discard tasks). 


3. A large percent of the items must survive to the 
wearout age or life-limit. 


4. A safe life-limit for an item must be established 


at an age below which relatively few failures are 
expected to occur. [Ref.10] 
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MIL-STD-2173(AS) points out two key elements for time-based 
replacement, namely an increasing failure rate, i.e., (2) 
above; and a large rate of survival to the wearout age or 
life-limit, (3) above. An increase in the occurrence of 
repairable failures with system age will often tend to result 
in a corresponding increase in maintenance costs. A decrease 
in system reliability and availability will also occur. From 
these two viewpoints, optimal replacement intervals should 
exist to minimize costs and/or maintain the system above 
certain minimum reliability and availability requirements. 

A large percent of items must survive to the point of 
wearout, i.e., point of increasing failure rate, or life- 
limit, to make planned replacements an effective maintenance 
tool. Since a scheduled replacement is deemed to be more 
desirable than an unscheduled one, the opportunity to make 
replacements should be utilized fully. Of course, the 
resources associated with the logistics of planning a 
maintenance evolution may be wasted if a large percentage of 
the replacements are premature and occur prior to the planned 
interval. 

Several additional factors must also be considered in 
determining the applicability of time-based replacement. Even 


a system with a constant failure rate may warrant time-based 
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replacement if an increase in operating and maintenance costs 
occurs as the item ages, or if repairs at sea are less easily 
made. Failures may occur no more frequently with system age, 
but the nature and type of repairable failures may result in 
an increase in the cost of parts and labor as the system ages. 
Another consideration in a time-based replacement scheme is 
the scope of work required to replace the component and the 
maintenance requirements of neighboring systems. Some 
replacements may require extensive interference removal, 
elaborate pre-established plant conditions, and extensive 
post-installation testing. In such cases, common sense 
dictates combining maintenance actions requiring the same or 
similar conditions. Such scheduling may not coincide with the 
optimal interval to minimize the operating cost for every item 
or subsystem, but even if some compromise is required the 
overall savings could be substantial. This thesis does not 
address the problem of coincident replacement of sets of 
different subsystems that have age-dependent failure 


properties. 


D. APPLICABILITY OF THE CAM-DRIVEN RECIPROCATING PUMP 
The cam-driven reciprocating pump is a candidate for time - 


based replacement. The data analysis of Chapter II suggests 
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an increasing occurrence of repairable failures as the pump 
ages, thus meeting criteria (b) for applicability. As stated 
earlier, the cam-driven reciprocating pump is currently on a 
time-based replacement schedule of 36 and 60 months for 
Trident and non-Trident assets, respectively. The remainder 
of this chapter will examine available data for evidence that 


the data supports the current replacement intervals. 


E. DERIVING A COST MINIMIZING FUNCTION 

As stated earlier, a stochastic model of the pump's failure 
characteristics allows the use of mathematical methodology to 
predict and quantify future behavior. One such tool is the 
Renewal Reward Process. 

Recall that the Nonhomogeneous Poisson Process is used to 
model individual pump failure times. The NHPP is not a 
renewal process, but since all new pumps are assumed to be 
similar, the number of pump replacements to occur in (o,t) 
does constitute a renewal process. Let M(t) represent the 
number of pump replacements occurring in a system up to and 
including time t, and let L,, n 2 0, represent the interval 
time between pump replacements. The cost associated with each 
renewal is denoted by R,. It is assumed that {L,} and {R,} are 


sequences of identically distributed random variables; a 
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generic L,, or R,, is denoted by L, or R. Then let R(t) be the 


sum of all system costs incurred by time t, so 


MO 


R® -¥ R,, (3.1) 
nd 


Further, denote the expected values for R, and L, as follows: 
E[R] = E[R,] and E[L] = E[L,]. Then the following proposition 
holds: 

Proposition 3.1. [Ref.11] 


If E[R]< ~ and E[L]< », then 


RO | EIR) 
t El] 


f+ & 


Let a cycle denote an individual pump's life. The proposition 
states that the long-run average cost equals the average cycle 
cost, i1.e., repair and replacement costs, divided by the 
average cycle length, i.e., pump life. More precisely, the 


long-run average system cost equals 


E[cost incurred during a pump 's life] (3.2) 
E[pump life] . 
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Pump life is a function of the replacement interval and the 
pump's mortality. Thus, a replacement interval for minimizing 
costs can be found by minimizing the long-run average cost, as 


given by Equation 3.2, over different replacement intervals. 


1. Expected Pump Life 

Let X be a random variable representing the age of the 
pump at replacement and let T be the designated maximum 
replacement age, in months. Further, define a cycle to be the 
interval from pump installation to pump replacement as defined 
by actual replacement or complete overhaul. Then the cycle 
length, denoted by L, equals X if the pump fails and requires 
replacement prior to the T; otherwise, L equals T if the pump 
life is at least as long as the scheduled replacement 


interval. Thus the cycle length can be summarized by 


Beye ge (3.3) 


It is assumed that the probability a failure is 
repairable is constant, i.e., independent of pump age and the 


failure number; that is, a failure is repairable with 
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constant probability p. The data have 140 failures, of which 
126 were repairable. This results in a point estimate for p 
of (126/140) = 0.9. The approximate normal 95% Confidence 
Interval for p is (0.85, 0.95). Likewise, define the 
probability that a failure is not repairable as q, q=1-p. 
Brown and Proschan used a similar assumption in their 
imperfect-repair model. The Brown and Proschan model assumes 
that the mode of repair was based solely on external 
conditions and not on the condition of the system at 
failure. [Ref .12] 

Table 3.1 is a contingency table showing the distribution 
of failures, conditioned on failure number. The assumption is 


that the proportion of non-repairable failures, q, remains 


constant over the number of failures. Let n,; denote the 
number of pumps having at least i failures, LoS 1p 2s oO 
Further, define Oo, as the number of observed replacements 


and 0;, as the number of repairs for pumps with at least i 
failures. Define the total number of replacements and repairs 
as R,, and R., respectively. The hypothesis of a constant 
proportion of non-repairable failures is tested by comparing 
the observed and expected values under the null hypothesis. 
Note that the observations for failures 5 and 6 are combined 


due to the small sample size. The test results in a value of 
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x? = 1.13 with 4 degrees of freedom for a p-value of 0.8 < p 
<0.9. The hypothesis of a constant is accepted at the 5% 
significance level. The small data set does not provide for 
the most accurate test, but it does give some indication of 
the goodness of fit. So the assumption of a constant 


probability for a failure being repairable is not unrealistic. 


TABLE 3.1:TEST FOR CONSTANT REPLACEMENT PROPORTION 


ee 
wmiacomes | s [> [2{2]2 | a 
Tamas foo foto [s | > | ane 
z 


CHI-SQUARED VALURE: x? = 1.13, DF = 4 
P-VALUE: 0.8 < p < 0.9 
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Now define h(x) as the probability density function for a 


pump's life, X. Then h(x) is 


ha) - + OAL, "Gg - e Meg A(x). (3.4) 


nd 
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Similarly, the probability that the pump does not fail before 


the scheduled replacement, P(X 2 T), is 


00 


P(X2T)- ¥ e sony, Ag OD (3.5) 
nt n: 


Using the above information it is possible to now derive an 
expression for the expected or mean life of a pump with a 
designated replacement interval of time T. In view of (3.3), 
the expected life, denoted by E[L], is 

T cs 
E[L] - fe h(x)dx + ft h(x)dz. (3.6) 
0 x. 
Inserting Equation 3.4 and further simplification results 


in the following expression for mean pump life: 
T 
E[L]- fp e ACW A(x)gde + T e AM, (3.7) 
0 
The above general formula can be specialized to account for 


any parametric form for A(t). Often the formula must be 


numerically evaluated; closed-form expressions may not exist. 


2. Expected Cycle Costs 
The cost incurred over the life of a pump consists of 
two components; the cost of a new pump; and repair costs 
incurred until the next replacement. Thus the life cycle cost 
of the pump is influenced by the number of repairable failures 


occurring over the pump's life. The expected life-cycle cost 


can be represented as 


Eflife-cycle cost] - COST | + COST :E[number of repairs 
new repair 


where both COST,.,, and COST,.,i, are expected or mean values. 
As stated earlier, the pump's life can terminate in one 
of two ways; the pump may have a non-repairable failure prior 
to scheduled replacement, or the pump will be replaced per the 
schedule at age T. Let N(t) represent the number of 
repairable failures occurring in the system up to and 
including time t. Likewise let E[N] represent the expected 


number of repairable failures occurring during pump life. 


Then 


E[N] - EINCQIX<T]-PUX<T) + E[N(1)|X2T]-P (27). (3.8) 


The analysis will examine each case individually. 
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a. Case X <T 
The expected number of failures given that the pump 


does not live to the scheduled replacement age T is 


= T < 
BINCO<T)-PU<T) - Yo nfe “ne AGP)" qdx. (3.9) 
nd ‘A n| 


This simplifies to 


T 
E[N(X)IX<T]-P(X<T) = gp fA@e “AGMD (x)dr. (3.10) 
0 


Integration by parts results in the following closed form 


expression: 


E[NCOIX<T)-PW<T) - A - pA(T)e™ - ms Bade (3.11) 


b. Case X > T 
The expected number of repairs given the pump 


survives to scheduled replacement is 


E[N(1)|X2T]-PX27) - x ne MAC)” (3.12) 
wd 


This expression reduces to 
E[N(1)|X2T] P27) - A(Type “™, (3.13) 


Finally combining Equations 3.11 and 3.13 the 
expression for the expected number of repairs over the pump's 


life is 
EIN] - 3 - pA(Tye %™ - Sa AM. A(T)pe MO, (3.14) 


The mathematical expression for the expected cost 


incurred during the pump's life is 


BIR]- COST ,,- COST yl? - pA(Tye M™ - Pe-M™ . ACType M™ (3.15) 
q 


q 


As with Equation 3.8, the above general formula can be 


specialized to account for any parametric form for A(t). 


3. Long-run Average System Cost 
The long-run average system cost can now be expressed 
in terms of Equations 3.7 and 3.15. Let z(T) denote the long- 


run cost average for the replacement interval of length T. 


Then 
COST, 44* COST pa,» | 2 - pA(T)e MP - Po, A(Type A 
2(T) = = 
f eG Aixjgde + T er 
° (3.16) 
Inserting the following parametric expressions for A(t) 
and A(t): 
e a 
A(x) 7a Bey) A(x) -0 BF, 

and p = 0.9, the optimal replacement interval is found by 


minimizing z(T), displayed in Equation 3.16, over T. 


F. DETERMINING THE OPTIMAL PUMP REPLACEMENT INTERVAL 
The first step is to determine representative values for 
the average replacement and repair costs. Here the 


replacement cost is simply the cost of a new pump which is 
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approximately $111,000, excluding the associated installation 
costs. The repair cost is a function of failure type. The 
failure type, however, may be affected by pump age. Some 
expensive-to-repair failure types may dominate later on in the 
pump's life, which can have a large effect on the average 
repair costs. Figure 3.1 exhibits the sample distribution of 
pump failures by failure type as the pump ages. Note, 
however, that the failure data is not adjusted for the number 
of operating pumps. Further analysis should include an 


estimation of failure rate for the specific failure types. 


DISTRIBUTION OF FAILURE TYPE 
(OVER PUMP LIFE) 
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Figure 3.1: Distribution of Failures by Type 
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This leads to a more useful representation of repairs 
cost: a moving average over pump age. Using information on 
material costs for repairs from the 3-M System, a rough 
estimate of repair costs for the different repair types was 
found. The 3-M data is sketchy, however, and information was 
not available for all repair types. 

The construction of the moving average involved ordering 
the failure times and assigning and an average repair cost to 
each based on the repair type. The moving average consisted 
of the average repair costs for ten consecutive failures. 
Figure 3.2 show this average cost as a function of time. Due 
to the limited information on repair costs, no attempt was 
made to further quantify the relationship of repair costs with 
age. No obvious trend is evident from the moving average, so 
for a conservative estimate of repair costs, the final average 
occurring at month 61, $4900, will be used as a basis for 
replacement interval determination. The repair cost estimate 
deals only with material costs. Other factors not captured in 
this cost estimate include labor costs and overhead costs. 
The 3-M System contains limited data on labor manhours for 
repairs, but again this information is sketchy. Similarly to 
the material cost projections, a ten point moving average 


was constructed to assess the changing maintenance 
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REPAIR ACTION MATERIAL COST 


(10 FAILURE MOVING AVERAGE) 


MATERIAL CO8T (8) 


17.2 19 23 25.2 32.3 36.2 
PUMP AGE (MONTHS) 


Figure 3.2: Repair Cost Moving Average (10 pt. Avg) 
requirements as the pump ages. No attempt was made to assign 
a dollar figure to these labor projections due to the 
complexity of the task and the lack of quality data. Future 
research, however, should be directed at developing such a 
cost relationship to be combined with material costs for a 
comprehensive evaluation of trends in Maintenance costs. 
Figure 3.3 show the resulting graph of the moving average 
for expended manhours for repairs as the pump ages. The graph 


shows a possible increase in labor requirements with pump age. 
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As stated earlier, no attempt is made to incorporate this data 


into the cost estimation. 


REPAIR ACTION LABOR REQUIREMENTS 
(10 FAILURE MOVING AVERAGE) 


Figure 3.3: Moving Average Repair Labor Requirements 
(10 pt. AVG) 


Determining the replacement interval is now simply a 
Matter of minimizing z(T) over T. A Mathcad 3.1 program was 
used to evaluate and graph z2(T) over a range of T. The 
optimal replacement interval is the minimum point on the 
graph. The optimal value of T is rounded to the nearest 
month. Table 3.2 contains the results for the MLE of oa and 


8 as well as for their 95% Confidence Bounds. Figure 3.4 is 


a graph of z(T) for the MLE of a and 6. Calculations are made 
for both the final moving average value of $4900 and the 
overall average repair cost value of $3000. Calculations are 
likewise made using the 95% Confidence bound derived in 
Chapter II. 


TABLE 3.2: OPTIMAL REPLACEMENT INTERVAL 


Coeu. = $111,000 REPLACEMENT INTERVAL 
= 0.9 (MONTHS ) 


a Ey 2 


foe) 
ies) 


“I 
“I 


The results indicate an optimal replacement interval of 
111 and 128 months for the average repair costs of $4900 and 
$3000, respectively. These figures as well as the calculated 
interval using the confidence bounds are well above the 
current replacement intervals of 36 and 60 months forTrident 


and non-Trident platforms, respectively. The failure and cost 
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data, therefore, may not justify the current Trident and non- 
Trident platforms, respectively. The failure and cost data, 
therefore, may not justify the current replacement interval on 


the basis of minimizing maintenance cost alone. 


LONG-RUN AVERAGE COST 


LONG-RUN AVG COST (z(t)) 


96 104 112 
REPLACEMENT INTERVAL (MONTHS) 


Figure 3.4: Long-Run Average Cost 


Recall, however, that only material costs are used for the 
repair cost estimate. Other costs associated with repair 
include the cost of labor and shipyard facilities. Even so, 
an average repair cost of approximately $32,000 would be 


required for an optimal replacement interval of 60 months 


using the MLE parameters. The lowest average repair cost 
required for a 60 month interval is $11,000 and occurs at the 
bounds of a = -2.87 and ® = 0.033. Other considerations not 
accounted for here, however, may be included in the current 
replacement policy. 

In developing a replacement policy for the cam-driven 
reciprocating pump, one must not lose sight that the pump is 
installed in a warfighting ship, namely a submarine. Cost is 
not and should not be the lone factor in establishing 
maintenance policy. System reliability and its effect on 
overall mission accomplishment must be a strong consideration 
in any decision. To aid decision makers, the following 
estimates of pump performance are calculated using the model: 
expected number of failures for a specific pump age, the 
expected failure times, and pump reliability for a specified 


mission duration. 


Since a NHPP with rate A(t) is used to model the failure 
characteristics of the pump the mean value function, Equation 
1.3, is used to calculated the expected number of failures 
over pump age. Figure 3.5 shows the expected total number of 
failures as the pump ages if the pump is never replaced. 


Recall that each failure is assumed to have a probability of 
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q = 0.10 that the failure is non-repairable. The sequence of 
failures, thus constitute a geometric distribution. The mean 
number of failures before a required replacement is (1/q) or 
10 failures, i.e., a replacement will be required for the 


tenth failure. 


EXPECTED NUMBER OF PUMP FAILURES 
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= 
a 
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Figure 3.5: Expected Number of Pump Failures 


The average times to failure can be calculated using a 
simple simulation program. Since the NHPP possesses a 
continuous mean value function A(t), Ginlar (1975) discussed 


the following recursive algorithm for generating a sequence of 


arrival times, t,,t.,...t,: 


1. Calculate the expectation function A(t) 
and its inverse A‘(t): 
t t 


A(t) = [P@dae s fe oBe ay 
0 0 


X- A(p - aC Bry) 
XB 


e” 


ACD - am 1) 


2. Generate a random variable U ~ U(0,1). 
3. Set t', =t'y.. - In(U). 
4. t, = A“[t',].[Ref.13] 

Table 3.3 contains the mean pump age at failure using 500 
replications of the above algorithm for simulated failure 
times. 

The decision to replace or repair a pump may be guided by 
certain minimum requirements for mission reliability. Recall 
from Chapter I that the reliability of the system at any time 
t, depends on the age at which the most recent failure 
occurred. Equation 1.5 gives the expression for the 


reliability of the system at that time. Figure 3.6 is a 


graph of the probability that a pump having a repairable 
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SIMULATED PUMP FAILURE TIMES 


TABLE 3.3: 


BETWEEN 
FAILURES 


failure at age t will successfully complete a mission of 
three months duration without experiencing a failure. Figure 
3.7 shows the expected number of failures that will occur over 
a 3 month mission given the pump is age t at mission 
commencement. This is computed using Equation 1.3 for pump 
ages (t,t+3). 

Both Figures 3.6 and 3.7 indicate that the replacement 
interval of 111 months, based entirely on material costs, may 
not be desirable due to the poor performance of the pump in 
terms of mission survivability, and expected pump failures 
during a mission. The current replacement interval of 60 
months may be acceptable based upon minimum mission 
reliability standards. No attempt here is made to define 
those standards. The replacement interval of 36 months is 
probably premature and should be extended at least to 60 
months to coincide with the policy for non-Trident submarines. 

Table 3.4 provides the estimated mission survivability and 
expected number of failures for a three month mission for 
pumps with mission completion ages of 36, 60, and 111 months. 
Also listed is the long-run average cost for a time-based 
replacement schedule with the associated pump ages. This 
table illustrates the trade-offs in cost and reliability that 


must be resolved for an effective maintenance policy. All 
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RELIABILITY FOR A 90 DAY MISSION 


PROBABILITY OF WO FAILURES FOR 90 DAYS 


4 
POMP AGK AT LAST FAILURE (MONTHS) 


Figure 3.6: Reliability for a 90 Day Mission 


EXPECTED WOMBER OF FAILURES DURING 90 DAY MISSION 
FOR PUMP AGE 7 


PUES OF EEPECVED PUD FAI LEOeS 


Figure 3.7: Expected Number of Failures During a 90 
Day Mission as a Function of Pump Age 


decision makers should understand the possible consequences 
of the proposed decisions. The use of probabilistic modeling 
is one way to evaluate the trade-offs and consequences of 
different maintenance policy. 


TABLE 3.4: MISSION SURVIVABILITY AND EXPECTED 
FAILURES FOR A 3 MONTH MISSION 


PUMP AGE PROBABILITY | EXPECTED LONG - RUN 
AT MISSION OF NUMBER AVERAGE 
COMPLETION NO PUMP OF PUMP COST 

(MONTHS ) FAILURES FAILURES ($ /MONTH) 


The calculations used in the preceding discussion are 
relatively simple to perform. Many similar calculations can 
also be made to address specific reliability requirements and 
different measures of effectiveness. The goal of this 
analysis has been to provide decision makers with a 
comprehensive assessment the cam-driven reciprocating pump's 
performance. Such information is necessary in the formulation 
of maintenance policy designed to address both economic and 


reliability concerns. 


IV. MODEL ASSESSMENT 


A. VALIDITY OF THE DATA 

The formulation of a specific NHPP to model the failure 
characteristic of the cam-driven reciprocating pump was 
accomplished with a MLE calculation using actual failure data; 
data obtained primarily through the Navy's 3-M system. The 
accuracy of the model, therefore, depends largely upon the 
degree to which the data accurately represents actual system 
performance. 

The current data collection method is far from perfect. 
The 3-M system itself suffers from many flaws. Part of the 
problem is inherent in the 3-M system, itself, and part is due 
to the fleet's attitude toward the system. When a failure 
occurs, the reporting process involve a crew member, normally 
junior enlisted personnel, filling out an OPNAV Form 4790/K. 
Some of the information to be included on the 4790/K are the 
equipment identification code, the date of failure, symptoms 
of failure, cause of failure, required repair parts, and 


required repair hours. 


The 3-M system is designed to track failure data for 
everything from mechanical and electrical systems to fuel 
storage tanks. In doing so the system is so large and generic 
in reporting criteria that problems arise from a lack of 
standardization in recording specific failure information for 
individual pieces of equipment. This is often exemplified in 
general and nonspecific entries for failure symptoms and 
causes which can lead to confusion in reconstructing the 
actual equipment performance. The data can be further 
confounded by improperly entered identification codes for 
equipment and repair parts as well as incomplete entries. The 
3-M system also does not capture all work performed by 
shipyard personnel during non-availability periods. All of 
these factors act to cloud the picture of true system 
performance and thus reduce the accuracy of any data analysis. 

The data used in the formulation of this model may not be 
totally accurate in its portrayal of the cam-driven 
reciprocating pump's maintenance history, but it is currently 


the only viable source of data available for analysis. 


B. VALIDITY OF THE ASSUMPTIONS 
In formulating the model, eight assumptions are made as to 


the characteristics of the cam-driven reciprocating pump 
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regarding failures, repairs and operation of the pump. Some 
of these assumptions are easy to accept; others require some 
discussion. First, this analysis assumes all pumps regardless 
of the submarine in which they are installed, experience 
roughly the same operating cycle. Pump operation during a 
submarine deployment will be similar between individual 
vessels. The deployment schedule, however, will differ 
between submarine platform and individual units. The long-run 
average deployment time is assumed to be the approximately the 
same for all submarines in the study. This, therefore, is a 
reasonable assumption. 

The assumptions of independence between consecutive 
failures and repairs returning equipment to full operation are 
related. In reality these assumption are not always true. 
One failure can cause a subsequent failure at a later age. By 
the same token, the act of affixing repairs has been known to 
cause future failures. These failures may be totally 
unrelated to the previous failure, but were caused by 
improperly restoring the equipment. Likewise, a repair may 
not fully repair the problem and the pump is left ina 
condition below full operating capacity. This same argument 
can be made for a complete overhaul of the pump in that the 


overhauled pump is not as 'good-as-new'. Postrepair tests and 
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procedures, hopefully, identify and correct such faulty 
repairs. In regards to this data set, suspect failures 
following a repair were carefully scrutinized to catch such 
double-failures. In the end, however, all such dependence 
between failures and incomplete repairs cannot be sifted from 
the data set. Therefore, they will cause some loss of 


accuracy to the model. 


C. ACCURACY OF THE MODEL 


-beware of mathematicians and all those who make 
empty prophecies. The danger already exists that the 
mathematicians have made a covenant with the devil to 
darken the spirit and to confine man in the bonds of Hell. 

St. Augustine 


St. Augustine may not have been talking about modelling the 
failure characteristics of a cam-driven reciprocating pump, 
but he does provide some wisdom for modelling in general. In 
fitting a mathematical model to characterize the failure 
behavior of any piece of equipment, it is naive to think that 
the model can perfectly predict the future performance. This 
would be true regardless of the system or the quality of the 
data. The results of this thesis, therefore, must be regarded 
in this light. This does not imply that any such model has no 
merit, but common sense and good engineering principles must 
be incorporated with any analytical results. The 
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incorporation of such information with the modelling results, 
provides a sound basis for making policy decisions. 

The previous sections discussed several areas that may 
contribute to providing inaccuracies in the model, namely with 
the accuracy of the data and the validity of the assumptions. 
Another factor that must be noted is that the original data 
set only provides information on pump performance up to about 
61 months. The majority of the observations actually occurs 
below 40 months. The model, however, extrapolates this 
information into performance predictions far and above this 
age. Any model inaccuracies will be magnified in this region. 
Also the use of this particular NHPP is only one parametric 
estimate of performance, many other more accurate models may 
exist. 

Given all the model inaccuracies, does this model provide 
any information to policy makers? Yes! The formulation of 
this model and the associated data analysis have revealed 
several important features of the pump 's failure 
characteristics. First, the pump exhibits wearout as noted by 
an increasing failure rate with age. Secondly, a time-based 
replacement decision of 60 months cannot be based solely on 
economic reasons. Thirdly, the model provides an estimate of 


the expected failures, failure times, and mission reliability 
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as the pump ages. This information, while not being totally 


precise, is pertinent to any decision maker evaluating 


maintenance policy. 
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V. CONCLUSION AND RECOMMENDATIONS 


As stated in the Introduction, the goal of this thesis has 
been to address the analytical needs of the Submarine Force 
necessary for a thorough assessment of current maintenance 
policies. Towards this end, an analysis of the failure 
characteristics of a cam-driven reciprocating pump has been 
conducted to demonstrate the degree of analysis required for 
a comprehensive assessment of equipment performance. Such 
analysis is necessary to assist decision makers in formulating 
maintenance policy especially when the Navy is faced with the 
reality of a shrinking budget. Decisions, however, must not 
be made strictly based on monetary measures. A thorough 
understanding of the consequences of any decision in terms of 
system performance and reliability is essential. 

This analysis has been an extension of the work preformed 
by Mr. Richard Youngk at SMMSO in evaluating the operation of 
the cam-driven reciprocating pump. His assistance and the 
support of SMMSO was vital in the completion of this work. 

The results of this thesis are not perfect. This analysis 


used one particular model, namely a NHPP with a specific rate 


function, to model failure characteristics. Other 
probabilistic models exist and a different model could 
possibly have provided a more accurate representation of pump 
performance. It could be advantageous to compare the results 
of this thesis with some of the other models. Additionally, 
one slight modification to the current model is recommended; 
perform a logarithmic transform on the failure times and 
recalculated the MLE. This would provide a smaller increase 
in the failure rate over time and may improve the model's fit. 

As stated earlier, the accuracy of this model, and in 
general any model, is to a large degree reflected by the 
quality and quantity of data. The modelling process is not 
complete; as more information becomes available, the model 
should be updated. Likewise, future work should be devoted to 
improving the data collection system. If the 3-M system is to 
provide quality information, improvements must = occur, 
specifically in the timeliness and accuracy of reporting 
performance data. Many independent data tracking programs 
have grown out of frustration with the current system. Such 
systems while providing quality information, place additional 
burdens on analyst and fleet personnel. The 3-M system needs 
to re-examine its purpose and assess how well it addresses the 


analytical needs of the Navy. 
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The evaluation of the cam-driven reciprocating pump should 
continue. New, quality data will improve the analysis and 
better quantify pump performance. Specifically, comprehensive 
cost and labor estimates will more clearly define the changing 
Maintenance requirements of the pump as it ages. 

This analysis has looked at only one small aspect of pump 
performance. Many additional areas merit research. One area 
centers around identifying differences in pump performance 
between platforms and individual submarines. The recognition 
and investigation of such differences could identify specific 
operating and maintenance practices unique to certain 
submarine that either enhance or degrade system performance. 

Another area for further study involves analysis of the 
failure patterns by specific failure type. The thesis model 
considered failures as one entity when in fact many different 
failure types exist. If a pattern follows the occurrence of 
certain failures, then future failures and failure types could 
possibly be predicted based on previous failure. This would 
aid in planning maintenance, but more importantly it would 
steer investigation toward why such tendencies exist which 


could lead to equipment or procedural modifications. 


The concept of failure prediction and prevention through 
monitoring also deserves attention. This again involves an 
analysis of failures by failure type, but additionally 
requires the identification of any precursors or indicators of 
impending failure. Such monitoring and predictive analysis 
already exist for many system. Further research could explore 
the value of certain precursors including determination of the 
expected time to failure given a specific precursors exist. 
This could aid in evaluating the effectiveness of monitoring 
procedures and possibly provide an optimal interval for 
planned monitoring. 

One final area of research is that of time-based 
replacements for multiple systems. This concept was discussed 
briefly in Chapter III in which maintenance for different 
systems may be scheduled to coincide if significant 
interference removal or abnormal plant conditions are required 
for both. This thesis only addressed one system, but multiple 
systems could be analyzed similarly to determine the most 
effective replacement schedule to maintain the readiness of 
all systems involved. 

This thesis has demonstrated the benefits of quantitative 
analysis and stochastic modelling in evaluating equipment 


performance. Such analysis are necessary to make the tough 
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decisions on maintenance policy in the wake of a shrinking 
budget. SMMSO, and specifically Mr. Richard Youngk, has taken 
on the task of incorporating these and similar concepts into 
current maintenance planning. It is hoped that policy makers 
will use these mathematical tool to make the well informed 
decisions necessary to maintain the U.S. Navy in top materiel 


readiness. 
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APPENDIX A. PUMP FAILURE DATA 
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FINAL EXCESS TOTAL 
STATUS LIFE LIFE 


Te 
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43 
41 
21 
10 
32 
13 
20 

9 
27 
22 
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43 
41 
21 
10 
32 
13 
20 
10 
29 
24 
26 
14 
13 
27 
10 
10 
11 
14 
13 
23 
17 
25 
33 
27 
36 
21 
22 
14 
52 
40 
37 
28 
53 


FINAL EXCESS TOTAL 


PUMP T, 7, T, T, T, T, STATUS LIFE LIFE 
34 14 34 Cc 16 50 
35 16 38 R 0 38 
36 18 41 Cc 18 59 
37 26 61 Cc 4 65 
38 11 12 15 Cc 11 26 
39 4 6 21 Cc 15 36 
40 5 20 23 Cc 30 53 
41 10 17 26 R 0 26 
42 17 32 33 Cc 3 36 
43 36 37 41 Cc 3 44 
44 16 23 43 R 0 43 
45 17 25 43 Cc 1 44 
46 18 24 56 Cc 7 63 
47 2 8 9 10 R 0 10 
48 15 17 22 25 Cc 1 26 
49 15 32 38 39 Cc 17 56 
50 27 32 37 42 R 0 42 
51 3 27 42 48 Cc 2 50 
52 4 19 52 52 Cc 8 60 
53 6 15 17 20 24 Cc 0 24 
54 9 12 16 26 26 Cc 3 29 
55 4 6 14 25 34 Cc 3 37 
56 5 16 18 36 43 Cc 2 45 
57 4 10 25 28 44 R 0 44 
58 19 31 42 46 49 Cc 4 53 
59 12 36 43 44 51 Cc 11 62 
60 11 22 25 34 35 37 R 0 37 
61 13 17 30 32 36 40 Cc 13 53 


LEGEND : 

T= FAILURE TIME, IN MONTHS, FOR FAILURE NUMBER i 
C INDICATES PUMP LAST OBSERVED OPERATING 

R INDICATES PUMP WAS REPLACED AT LAST OBSERVATION 


NOTES : 

ALL TIMES ARE ROUNDED TO THE NEAREST MONTH 
TOTAL FAILURES: 140, TOTAL REPLACEMENTS: 14 
TOTAL EXPOSURE: 2020 
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APPENDIX B. DERIVATION OF THE MAXIMUM LIKELIHOOD ESTIMATOR 


k = index for the number of pumps in the sample set 

j = index for the individual pumps, j = 1,2,..k 

i = index for the failure number for pump j, i = 1,2,..N, 

n, = the number of failures for pump j 

f, = time of the last observation for pump j 

*Note: time is equivalent to pump age in this analysis* 

t,, = time of failure i for pump j 

p(t,;) = probability that the failure occurring at time 
t,;,, is repairable. 

q(t,,;) = probability that the failure occurring at time 
t,, is not repairable and therefore requires 
pump replacement. q(t;,;)= 1-p(ti,;). 

F, = indicator variable indicating status of the last 

observation 
F, = 1 for non-repairable failure 
F = 0 for pump last observed operating 

A(ty,;) = rate of occurrence of pump failures for pump age 
t, .. The failure may or may not be repairable. 

A(t;,;) = the integrated value of the failure rate from 
pump installation to age t,,,; the expected number 
of failures(repairable) up to age tj; 

6 = vector of parameters, 6 = (a, 6) 

L(@) = Combined Likelihood Function for pump failures 

L,(@) = Conditional Likelihood Function for pump failures 


NOTE: If the last observation is a non-repairable failure then 


f. =t and F. = 1. 


| nj,3 
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Consider the likelihood function for one pump, i.e., the j*: 


L{0)- “PA, ptt, de EMOPMODN 6, ype). 


aa IG, | : (B1) 
Phy) 


Note that in (B1) the dependency of A, A and p on parameters 
6 is suppressed. 

For the case with the last observation being a non- 
repairable failure, f; = t,j;,; and F; = 1 so expression (B1) 


reduces to 


L{0)- e “2, drt, De EMP (6, YP (ts) 


(B2) 
0G, alta) 
Equation (B2) simplifies to: 
as Ue)" 
{Oy AGE MEDI) Meda} oe 
(B3) 


ne TT 10,306 {24)" 
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Combining the observations from all pumps: 


ieee “TTT La, ah 2 Hla.) (B4) 


fl él is 


To obtain the conditional probability density function, or 
pdf, of the observations given n,, j=1,2,...k, failures, 


divide L(8) by 


This results in the conditional likelihood function: 


Il [be If», ! 


L.@} jel id (B5) 


Ia)" 
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SPECIAL CASE: MODEL FOR INCREASING HAZARD 
In the remainder of this thesis we consider the following 


specific model, developed by Cox and Lewis (1966), pp. 45-54: 


AG)-e" 4 AG)» eh. 


Then the combined likelihood function, L(8), becomes: 


BS 8S, Leetole xt, 
4) «ele UY IB 5 Biter . (B6) 


ptt, 5 


Similarly the conditional likelihood function, L,(8), is 


at “TT ; 


L (8)=§ ——_2>__—_ (B7) 
fle’ J-}) "i 
Here 0 = (a, ®) and any parameters used to specify p(t;,;). 
Taking the logarithm of L,(8) results in 
kon k k k 
log L@} BY Y 4, :Yolog(y,)) log BY”, - Yon Joste™-1) (B8) 
fl id fl ji ja 
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To find the Maximum Likelihood Estimator for B, take the 
derivative of log L,(@), equation (B8), with respect to B, 
set it equal to zero, and solve for B. The resulting 


derivative is 


dog LAO) Sh ifn kos, 
—_—_ ifs as . 


ny - 
4B A ‘tl Bot 7 et ey 


(B9) 


Graphing Equation (B9) over a range of values for B and 
visually locating the zero intercept produces a rough 
estimate of the maximizing value for 6. This initial 
approximation can be further refined by a using a numerical 
method, such as Newton's Method, to solve for the root. This 
reduces the mathematical complexity by eliminating the need 
for deriving a closed form solution for the root if a closed 
form solution does indeed exists. 

A similar method is used to find the MLE of a. The 
Conditional Likelihood Function, L,.(8), however, cannot be 
used to solve for a as it was for 8 since the conditioning 
process removes @ as a parameter. Therefore, L(@), Equation 
(B6) must be used to find the maximizing value of a. Taking 


the logarithm of L(8), differentiating with respect to q, 
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and setting the equation equal to zero results in the 


following expression: 


flog LOY) ee hay. 
a ae aye )- 0, (B10) 


Solving the equation in terms of a results in the following 


expression: 


& = nj—#2_| (B11) 


Inserting the maximizing value of B into Equation (B11) 


produces the resulting maximizing value of a. 


CONFIDENCE BOUNDS FOR « AND £6 

The confidence bounds for the estimation of a and ® are 
derived using the Fisher information matrix to obtain the 
asymptotic variances and covariances of the MLE's. The 
Fisher information matrix is the composed of the negative 
second partial derivatives of the sample log likelihood. 


Inverting the Fisher information matrix evaluated at the MLE 
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for a and B produces an estimate of the variance/covariance 
matrix. The variance can then be used to derive the 
confidence bounds for the MLE's. 

The partial derivatives are derived using the log of 


equation (B6é) and are as follows: 


PlogL(0) - ee an 
eee). -1) 
5a? 8 » (e (B12) 
alogL.()_ atogL(8) - ey (e Pty : ety (f, By | (B13) 
da dp ap ac p? yl B yl f 
Plogh (8) = 2e “wr (Bh yy , 26“ (rp em 2 BF 
Bol Le as idl Me ce CL) ote by eS ty, 
a5 ms Yt ar ¥ Uf ee ¥ Ge (B14) 
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APPENDIX C. MATHCAD 3.1 PROGRAMS 


RII RAIA IIIA AIA II III AI IIASA IAI ASIII IIASA IASI ISIIS IAI IAS SAS AL ISI SISA IA ISOS ISIS IASI SAS IISA ISIS ISAI IIIA IS IAI ASIII IIR, 


MATHCAD 3.1 PROGRAM TO COMPUTE MLE OF BETA 


AKI IIIA IIA IOI II IAAI IAAI IAI INIA IA IAI AAI ISIS IIIA IIS SISI IASI ASI IIIS ASIN ISIS I AISI SSIS IIA III IISA SSSI SASSI SIN IIA 


j-=0.60 k -1.. 150 
n =READ(NU) tt, = READ(TTU) f, = READ(TFU) p, =—— 


1 (Nore de Ly ia ma) pe 


THE DERIVATIVE OF THE CONDITIONAL LOG LIKELIHOOD 
PLOTTED OVER VALUES OF BETA 


fI(B) 
B 
NEWTON'S METHOD 
N = 60 i=0..N aah 
y, 7-023 fil (Yo) =-7.88315 re 
es 0.0225818402 
y. 
ae a )) | - err, ly, - rs y =| 0.0225823081 
(¥) 0.0225823081 
0 
2 -=last(y) - 1 
n2 =3 
fl (Yaa) =0 
Y, = 0.02258 


80 


KRRHK KREIS IKI AERA IARI EER IARI IIIS IIASA IIASA III IISA AAAI ITAA AREER 


MATHCAD 3.1 PROGRAM TO DERIVE THE 95% Cl FOR MLE'S ALPHA AND BETA 


HRKKIIRIKIHAI IIR ARIA AIA AIA II ARIA AR IARITASRIIAI IIA IIIA IAI AIA IAAI AIA IASI RATA IIIS AAI ARE RARER 


j -0.. 60 


n, -= READ(NU) tt, -= READCTTU) f ‘= READ(TFU) 


Second partial derivatives of the unconditioned MLF 


daa(a,b)=-£-5) (53) aab(a,b) “tye 


J j j 
22 fprle ell ALI] ee gr 
dbb(a,b) := : y. “ai ae if ) ae pe (. rE: oe (. ).(6) 
j j j j 
alpha ‘= - 3.189 
beta := .02258 
daa( alpha, beta) = -139.441 are - dab(alpha, beta) \ | 
5 i 

dbb( alpha, beta) =—1.257°10 -dab(alpha, beta) - dbb( alpha, beta) 
dab(alpha, beta) =-3.56°10° 

0.026 -7,335-10 ¢ | 

MAT = - 3 
-7.335°10 4 2.87310 > 


95% CONFIDENCE INTERVAL FOR ALPHA 


LLA = alpha~ 1.96 (MAT, ,) = 
LLA =~3.504 


95% CONFIDENCE INTERVAL FOR BETA 
LLB := beta — 1.96- (MAT, : 2 


LLB = 0.012 


ULA ‘= alpha+ 1.96-(MAT, 4) e 
ULA =-2.874 


ULB = beta + 1.96: (MAT, ,)° 


ULB = 0.033 


HERI I IAI IAI PARA AAAI AAIS APA AA ADI AAAI ASIANA IAAI ISSN SSSI ISD IIA ISAS AAAAIASSI IS IIIA ISSA III ISI III IR, 


MATHCAD 3.1 PROGRAM TO DETERMINE THE OPTIMAL REPLACEMENT INTERVAL 
Sect te tet tc RIAA AAI AA AAAI NIAAA SAAS SASS AIA S SANA AI AIA ISI IAID ASSIS AAI IASI PISS I SSR PASISSS A IASAAPISA ASAI SIAR 


Expected life of the pump (ELife): 


T = timed replacement interval in months 
p = probability that a failure is repairable 
q = probability that a failure is not repairable and requires pump replacement 


i =1..200 Tei 
126 
a'=-3.189 b:=0.02258 p =—— Shop 
P7140 ‘ 
e* /b b 
act) =©-(e*- 1) dA(t) =e" 


b 


ttf = case in which pump is replace before scheduled interval T 
ttT = scheduled replacement at time interval T 


T. 
‘ -A(T.\-q 
ttf. | q:xe 841 dcx) a tT, = Tye ( , f(x) 2 gx MO)" GA(x) 
0 

ELife. = ttf, + tT 


Expected rewardicost : 
CReplace ‘= 111000 


CRepair ‘= 4900 


-A(T.) - -A(T.)- 2 
ECost. = CReplace + CRepair: A(T,)-pre (") ibe A(T,)-pe oe Pie 
1 1 q 1 q 
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Reward / Renewal Process: Long Run Average Reward (z): 


3000 min, = 10000000 


Zy°= 1000000000000 


j =1..300 


i, voi (2, . (8) 


2798 


pid n2 ‘= last( min) 


lf 


Zz. . =2.22645726:10° 
n2-1 


aes 2, = 2.22641167°10° 


Zz. . =2.2264619*10° 
n2+1 


2192 n2=111 


FEES SIO III III IIS TTT IR III III IAI TI IIIS SISSON ISHII AIST SIS IS III I SAI I IT III IA IIT 


MATHCAD 3.1 PROGRAM TO SIMULATE PUMP FAILURE TIMES 


FEE IIIA III I III II III IT IIIA I SOSA SAIS ISIS IIIS ISI ISS IIIS ISIS ISS SSI ISS I ASI SISA ASIII IIIT III 


Simulation of failure times using method of Law & Kelton, 1982, pp. 509-510. 


ORIGIN=0 n:=500 = a:=-3.189 


i=1.14 


uj j--™d(s) 


jozlun 


eg Ranges he 


b :=.02258 


tpy | =O i=l 


fail9 = (7) — 


fail13 = a 


CALCULATION OF MEAN TIME TO FAILURE FOR FAILURE i 


i= {42 
faill: 


J 


avg, ce 


Dini, 


avg4 ed = 
yiai, 

a 

ave, J 2 


n 


avg, = 


aves 


aVBg : 


fail2 := (7) nas 

fail6 -= (:1) 2 
6 

fail6 «= «Qn : 


wie 


faill4 aa’ 


Disa 
; n 


disaits 
-j 
n 


dosails 
=_J = _ 


ave'=(0 16.668 29.288 39.841 
(88.276) (91.262) (94.209) 


Pree eae 
fail7 := (3) — 
Pers (a 


1 
faill] a . 


Disits 


J 


avg, := 
n 
yBils, 
- J 
avee = J h 
Vail, 
- J 
aVBo =J 
n 


Disilto, 


avB19 ~ 


Dail, 


aia = 
aite=() 
aie '=() 


faill2 age 


n n 


J 


n n 


isi 4 


n 


48.154 55.474 61.885 67.59 72.539 76.899 81.285 


84.829 


CALCULATION OF THE STANDARD ERROR FOR THE SAMPLE MEAN TIME TO FAILURE i 


jozlun 
r 5 r ae 
> (ail; = avg)) : > (tail2; 7 avg) 2 
samplestd | = he eaeeemeeae samplestd., =: eens 
n-1 n—1 
r ; 2) 2S : 5 
yy (dail; ~ avB,) »y (fail4, - avg,) 
samplestd., = ne samplestd , = j 
n—1 n—-1 
4.5 5 
by (Bails; = avgs) f YX (file; - avg¢) 2] 
samplestd, eo ene ee samplestd, := j 
n-1 n-1 
y 15 ; 2 ) 
(fail, ~ avg7) y (ail8, - avgg) 
samplestd, «= . ae samplestde -= j 
n- 1 n- | 
r 7.5 4.5 
> (ail9, - ave,) 2 by (tail10, — avg,,) ci 
samplestdy *=] J samplestd  , °=| 
a ae ged 
a, 5 r 4.5 
oD (Bill; ~ avg, ) 4 > (lail12; ~ avgy9) z 
samplestd,, = cares aan samplestd 5 |= a 
n— L n— 
r -) re] 
> (Bailt3; ~ avg,,) > (Bil14; ~ ava, 4) 2 
fo es ee lestd,, =| eee 
samplest 13 a samplest 14 I | 


stderror ‘las [sens 


(va) 


stderror?=(0 0.602 0.672 0.691 0.681 0.659 0.634 0.624 0.589 0.555 0.539 0.517 0.501 
(0.489) (0.479) 
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